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We determine the critical Zeeman energy splitting of a homogeneous Fermi superfluid at unitary 
in terms of the Fermi energy ep according to recent experimental results in LKB-Lhomond. Based 
on the universal equations of state for the superfluid and normal phases, we show that there exist 
two critical fields H c \ and H C 2, between which a superfluid- normal mixed phase is energetically 
favored. Universal formulae for the critical fields and the critical population imbalance P c are 
derived. We have found a universal relation between the critical fields and the critical imbalances: 
Hd = 7^ep and H C 2 = (f + *yP c ) 2 ^ 3 Hci where £ is the universal constant and 7 is the critical 
value of the chemical potential imbalance in the grand canonical ensemble. Since £, 7 and P c have 
been measured in the experiments, we can determine the critical Zeeman fields without the detail 
information of the equation of state for the polarized normal phase. Using the experimental data 
from LKB-Lhomond, we have found H c \ ~ 0.37ep and H C 2 — 0.44tF. Our result of the polarization 
P as a function of the Zeeman field H/ep is in good agreement with the data extracted from the 
experiments. We also give an estimation of the critical magnetic field for dilute neutron matter at 
which the matter gets spin polarized, assuming the properties of the dilute neutron matter are close 
to those of the unitary Fermi gas. 
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I. INTRODUCTION 

While Bardeen-Cooper-Schiffer (BCS) superconduc- 
tivity/superfluidity [f| in Fermi systems was investigted 
more than 50 years ago, the main scientific interest in 
recent experiments of two-component ultracold Fermi 
gas is to create a new type of Fermi superfluid in the 
BCS-Bose- Einstein condensation (BEC) crossover 
At the Feshbach resonance where the s-wave scattering 
length a s diverges^ novel type of Fermi superfluid has 
been observed p-Q- ^ n * ne dilute gas limit, where the 
finite-range effect of the interaction can be neglected, the 
only length scale in the many-body problem is the inter- 
particle distance. In this so-called unitary Fermi gas, 
any physical quantity can be expressed in terms of its 
value in the non-interacting case multiplied by a univer- 
sal constant [§]. For example, the energy density can 
be written as £ — ££q, where £ is a universal constant 
and £q the energy density of the non-interacting Fermi 
gas. A possible realization of such a unitary Fermi super- 
fluid is the dilute neutron matter [Tol— fl~2l| which exists in 
neutron star crust, since the neutron-neutron scattering 
length a nn ~ — 18.5fm is much larger than the effective 
range of the nuclear force and the inter-particle distance. 

In addition to the idealized case where fermion pairing 
occurs on a uniform Fermi surface, the effect of pure Zee- 
man energy splitting Ez — 2H (Here H — hbB with the 
effective "magnetic field" B and the magnetic moment 
/ie) between spin- up and -down electrons in the BCS su- 
perconductivity was known many years ago [j~3l - fl7l |. At 
a critical Zeeman field or the so-called Chandrasekhar- 
Clogston (CC) limit H c = 0.707A where A is the 
zero temperature gap, a first order phase transition 
from the gapped BCS state to the normal state oc- 



curs. Further theoretical studies showed that the in- 
homogeneous Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) 
state |16l If 7] may survive in a narrow window between 
H c and i?FFLO = 0.754A . However, since the thermo- 
dynamic critical field is much smaller than the CC limit 
due to strong orbit effect [l3| , it is hard to observe the CC 
limit and the FFLO state in ordinary superconductors. 

Recent experiments for strongly interacting ultracold 
Fermi gases give an alternative way to study the pure 
Zeeman effect on Fermi superfluidity [f8|, Ojj . The atom 
numbers of the two lowest hyperfine states of 6 Li, de- 
noted by Nf and N±, are adjusted to create a popula- 
tion imbalance, which simulates effectively the Zeeman 
field H in a superconductor. At the unitary point, phase 
separation phenomenon between the unpolarized super- 
fluid and the polarized normal gas, predicted by early 
theoretical works [20L [2l]| and Quantum Monte Carlo 
(QMC) calculations |22j |. has been observed [H, [Tg| . 
However, the evidences for the FFLO and the breached 
pairing states |23l |24| have not yet been found at uni- 
tary. The problem of unbalanced pairing in strongly in- 
teracting Fermi systems is also relevant to nuclear mat- 
ter [25] ] and quark matter (26j which may exist in compact 
stars. Since the neutrons carry a tiny magnetic moment 
fi n ~ 6.03 x 10- 18 MeV-G _1 (in natural units), the su- 
perfluid state in dilute neutron matter can be destroyed 
when the magnetic field reaches a critical value and the 
matter becomes spin-polarized 27]. Since the effective 
interaction in dilute neutron matter is quite close to uni- 
tary, the experimental data from cold atoms can help us 
determine the critical magnetic field for dilute neutron 
matter. 

Due to the universality at infinite scattering length, 
the CC limit of a unitary Fermi superfluid should be a 
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universal constant. For cold atom systems, the chemi- 
cal potentials for the j" and ], atoms are different. They 
can be denoted as [jtf and fi±. We define the averaged 
chemical potential fi = (/i-j- + /xj,)/2 and the Zeeman field 
H = (/i-j- — /i^)/2. For condensed matter systems, the 
Zeeman field H is induced by some external field. There 
exist two types of CC limit which can be directly mea- 
sured in experiments or QMC calculations corresponding 
to different ensembles : 

(l)Thc critical value 7 for the chemical potential imbal- 



ance 



(1) 



at which a first order phase transition between the su- 

pcrfluid and the normal phases takes place. 

(2)The critical value P c for the population imbalance 



P = 



N t -N± 



(2) 



above which the supcrfluid-normal mixed phase disap- 
pears. These two critical values have been successfully 
determined in experiments and QMC calculations in re- 
cent years |28l - t3l| . 

To make a closer connection with condensed matter 
systems and cold dilute neutron matter in neutron star 
crust, we should convert the results from cold atoms to 
usual variables of condensed matter systems, i.e., the crit- 
ical Zeeman energy splitting H c in terms of the Fermi en- 
ergy ep = (37r 2 n) 2//3 /(2M) (M is the fermion mass) for a 
homogeneous Fermi gas with fixed total density n. Obvi- 
ously, the ratio H c /ep should be a universal constant at 
unitary. However, from the simple mean-field approach 
in Ref. [33M35l j it was found that there should exist two 
critical fields H c \ and H C 2 in the BCS-BEC crossover and 
the supcrfluid-normal mixed phase is the energetically fa- 
vored ground state in the region of H c \ < H < H C 2- For 
H < H c i, the ground state is an unpolarized superfluid 
with balanced spin populations. The matter becomes a 
spin-polarized normal Fermi liquid for H > H C 2- The ex- 
istence of the two critical fields is somewhat like the type- 
II superconductors, but mechanism is different. Since the 
particle numbers iVf and N± are used as tunable parame- 
ters in QMC calculations and experiments, only the first 
order phase transition point 7 = (H/ fx) c and the critical 
population imbalance P c have been directly determined. 
The two critical fields H c \ and H C 2 in terms of the Fermi 
energy ep for homogeneous Fermi gases have not yet been 
determined by QMC calculations and experiments. 

In this paper, we will determine the two critical fields, 
H c \ and -ff C 2, in terms of the Fermi energy ep for a ho- 
mogeneous Fermi gas at infinite scattering length, based 
on the universal properties of the thermodynamics. In 
our study, we use the general equations of state for the 
superfluid and the normal phases [H, [13] ■ We neither 
assume that the normal phase is fully polarized [22L [38| 
nor adopt the non-interacting equation of state for the 



partially polarized normal phase in the naive mean-field 
approach [34| . Our results are therefore reliable once the 
QMC or experimental data are used as input. 



II. FORMALISM 
A. Equations of State and Phase Transition 

At unitary, we can construct the exact equations of 
state (EOS) in the grand canonical ensemble from the 
universality hypothesis. The pressure of the polarized 
normal phase (N) as a function of the averaged chemical 
potential fi and the Zeeman field H takes the form [9( 



H 



(3) 



where "Po(m) — |c/i 5 ^ 2 is the pressure of the non- 
interacting Fermi gas with c = (2Af) 3 / 2 /(37r 2 ), and Q{x) 
is a universal scaling function representing the strong 
coupling effect and can be determined from the exper- 
iments and QMC calculations. For the balanced case, 
we define G{0) = S,^ 3 ^ 2 and the universal constant £n 
for the normal phase should be larger than that for the 
superfluid phase, £n > representing the fact that the 
superfluid state is the ground state at H = 0. Experi- 
ments and QMC calculations have determined the value 
£ N ^ 0.51 - 0.56 [HEl. The total number density and 
the spin density imbalance, i.e., n^(/i, H) = dV^/dfJ. and 
m^(iJ,,H) = dV^/dH, are given by 



n N (n,H) = n {n) 



H 

2 / H 

m N (n,H) = n (fj)-Q' I — 
5 \n 



--Q' (- 

5 n V V- 



(4) 



where no(^) = c/1 3 / 2 is the density of a non-interacting 
Fermi gas. 

There exists a critical value 5q of the ratio above 
which the normal Fermi gas is in the fully polarized state 
(Npp) with a single spin component f, i.e., ttjn = «n 
or P = 1. In this case the Fermi gas should be non- 
interacting and the scaling function is given by Q{x) = 
+ For Pi j [i < Sq, the normal phase is partially 

polarized (Npp) with Tom < «n- While the mean-field 
theory predicts So — 1 [33j, it is 3.78 [HI from recent 
QMC calculations. Theoretically, the numerical value of 
<!>o can be determined by studying the impurity problem 
of a single J, atom immersed in a Fermi sea of t atoms [39l | . 

QMC studies [13, HH and experiments |29l - l3"i| have 
shown that the stable superfluid phase (SF) at unitary 
should be unpolarized. Hence the pressure does not de- 
pend on H explicitly and takes the well-known form 

Vsf(h,H)=V (v)C 3/2 - (5) 
The total density reads nsp(iJ, 7 H) — no(/i)£~ 3 / 2 . 
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FIG. 1: (a) Schematic plot of the equations of state for the 
normal and superfiuid states. The vertical dotted line shows 
the first order phase transition point H//i — 7. (b) The spin 
polarization P as a function of H/fJ,. As we expect from the 
first order phase transition, there is a jump from P — to 
P = P c at the first order phase transition point H/ fi — 7. 



With the QMC data £ = 0.42(1) and /? = 1.2(1), the 
authors in :22] obtained H c ~ 1.00(5) A !22j. However, 
in this case we found 7 = 1.22 < So — 3.78, using the 
same value of £ adopted in [22|]. This means that the 
normal phase should be partially polarized at the phase 
transition and hence the critical population imbalance P c 
should be smaller than unity. Therefore, the assumption 
of fully polarized normal state is not correct at unitary. 
In Fig. [TJ we demonstrate a schematic plot of the equa- 
tions of state for the superfiuid and normal phases, as- 
suming that the normal phase is partially polarized at 
the first order phase transition. 



B. Critical Zeeman Fields for Homogeneous 
System 

To make a closer connection with condensed matter 
systems and cold dilute neutron matter in neutron star 
crust, we should consider a system where the total den- 
sity n is fixed and the Zeeman splitting H is induced by 
a tunable external field. For neutron matter, the Zeeman 
splitting can be induced by a strong external magnetic 
field, since the neutrons carry a tiny magnetic moment. 
To determine the critical fields for homogeneous Fermi 
supcrfluids, we turn to the canonical ensemble with fixed 
total particle number density n = (2Mep) 3 / 2 /(37r 2 ). The 
Zeeman splitting H is treated as a real external field, and 
the conversion between particles in the states t and ^ is 
allowed, but the chemical potential \x is not a free pa- 
rameter. 

The chemical potential \i^(H) in the polarized normal 
phase is solved from the number equation 71n0uN) H) = 
n, and the energy density £n{H) = fiTst(H)n — 
Vn{hn{H), H) reads 



/'N 



5/2 



£0 (8) 



Considering the relation £n > £ and the fact that Q (x) 
is an increasing function of x, a first order phase tran- 
sition takes place when the pressures of the normal and 
superfiuid phases become equal, i.e., when H/[i reaches 
another universal critical value 7, which is determined by 
the following equation 

s(7) = r 3/2 . (6) 

For the case 7 > 5q, i.e., the normal phase at the phase 
transition is Npp , which was assumed in the early studies 
[H H, we have § (1 + 7 ) 5/2 = <T 3/2 - Combining with 
the relation for the excitation gap A(/i) = /3/j, at unitary 
for the balanced case (H = 0) with a universal constant 
P , we recover the CC limit derived in [22[ 

-^= 1 -( 2 ^-l) (7) 



with £0 — |ce^ 2 . At H — 0, we have /in(0) = Cn^f and 
£n(0) = £n£o- For nonzero H, the well-known relation 
£ = 3P/2 Q is broken, since the interacting energy with 
the external field H is included. In the fully polarized 
normal phase, we have /LIn(JJ) = 2 2 / 3 ep — H. The Npp- 
Npp transition occurs at H = 2 2 / 3 ep5o/(l + So). While 
the mean- field theory predicts Ho = 2" 1 / 3 e F ~ 0.794e F , 
we obtain Ho — 1.26e F from the QMC result So — 3.78 
& 

Solving the number equation nsp(fX,H) = n for the 
superfiuid phase, the chemical potential and energy den- 
sity are given by ^sf(H) = £e F and £sf(H) = ££ , 
respectively. At H = 0, there is the BCS instability 
£sf(0) < £n(0) which is numerically supported by the 
fact £n > £■ While £sf(H) keeps independent of H 7 
£n(H) should be a monotonously decreasing function. 
If there exists no heterogeneous mixed phase, a phase 
transition occurs at £^(H C ) = £sf- Thus H c can be 
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determined once G(x) is known. Furthermore, if we as- 
sume that the normal state at H c is fully polarized, i.e., 
H c > H , we obtain H c = §(2 2 / 3 -£)e F . Using the QMC 
data £ ~ 0.42, we get H c ~ 0.7ep which is smaller than 
Ho. According to the QMC data for the excitation gap 
Ao ~ 0.4 — 0.5eF pol . l4lj |. it is also in contradiction to 
the constraint H c < Ao which ensures that the superfluid 
phase is unpolarized. Therefore, the above assumption 
of phase transition between SF and Npp without SF-Npp 
mixed phase is not correct. 

The above simple analysis without heterogeneous SF- 
N mixed phase is not adequate since the first order phase 
transition should be associated with the phase separation 
phenomenon when the total density n is fixed. Now we 
take the heterogeneous SF-N mixed phase into account. 
In fact, due to the strong coupling effect in the BCS-BEC 
crossover, the critical field H c = H c (fi) = 7/i for the first 
order phase transition in the grand canonical ensemble 
splits into a lower critical field H c i and an upper one 
H C 2, and the SF-N mixed phase should appear in the in- 
terval H c i < H < H C 2. This can be understood by the 
fact that the chemical potentials for the normal and su- 
perfluid phases do not equate in the BCS-BEC crossover, 
which is different from the classical weak coupling case 
where the chemical potential is always set be equal to the 
Fermi energy e F . 

H c i and H C 2 can be determined by setting the chemical 
potential /i to be its value in the superfluid and the nor- 
mal phases, respectively. We have H c i = ^^sf{H c i) and 
H C 2 = 7Mn(#c2), where /xsf(-Hci) = Vn(H c2 ) can 
be obtained by the number equation ^ of the normal 
state at H — H c2 . We find 



Hn(H C 2) 



3/2 



70' (7)] 



2/3 ' 



(9) 



where we have used the fact £(7) = £~ 3 / 2 . Therefore, we 
arrive at the following model-independent expressions for 
the lower and upper critical Zeeman fields 



Hci = 7£ £ f, 



H, 



0* 



c-2 



[1 - ie/ 2 iG'(i)] 



2/3 



(10) 



The appearance of the mixed phase requires H c i < H c2l 
which gives rise to 



Q'il) > 0, 7 0'(7) < 2^ 3/2 - 



(11) 



Furthermore, since the superfluid phase is unpolarized, 
the condition H c \ < Ao should be satisfied, which leads 
to the fact 7 < /?. Therefore, the QMC result /3 ~ 1.2 
[ibl lilj gives an upper bound for 7. 

If we assume that the normal phase is fully polarized, 
the two critical fields can determined with only one pa- 
rameter £, 

/ 2 2/5 \ 



H C 2 — 



2 2/5 



2 4/ 
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£3/5 



■e F . 



(12) 



Substituting £ = 0.42, we obtain H& = 0.51e F and H C 2 = 
2A7ep. However, this is not the realistic case and is also 
not consistent with the QMC result for the pairing gap 
Ao = 0.4- 0.5e F S3, El. 



C. Properties of the Mixed Phase 

A SF-N mixed phase should appear in the interval 
H c i < H < H C 2, which essentially corresponds to the one 
observed in QMC calculations at < P < P c [28]. From 
the phase equilibrium condition Vsf{^tH) = "Pn(Mi H), 
in the mixed phase the ratio H/fj, keeps a constant 7 
and the chemical potential reads ^m(H) = H/j. The 
properties of the normal domain in the mixed phase de- 
pend on the value of 7. The QMC calculations [28| and 
experimental data support 7 < Sq, i.e., the normal do- 
main is partially polarized. The volume fractions of the 
superfluid and normal phases can be denoted by y and 
1 — y, respectively. y(H) is determined by the equation 
n = y(H)nsF(fJ-M, H)+[l- y(H)]n N (fj,M, H). Using the 
EOS for the phases SF and N, we find 



y(H) 



5£ 



-3/2 



27^(7) 



H 

H c i 



-3/2 



H C 2 

Hci 



-3/2' 



(13) 



The energy density of the mixed phase, £m(H) = jiyin - 
Vm (/im , H) , can be evaluated as 



£m{h) = \JL 

6 H c i 



2(H_ 
5 \Hc i 



3/2' 



t£o. (14) 



Since the normal domain possesses imbalanced spin pop- 
ulations, there is a nonzero global polarization P in 
the mixed phase, i.e., the system becomes spin-polarized 
when H > H c \. From the definition of P, we find 



H 

H c i 



3/2 



Hd < H < H c2 . (15) 



The mixed phase continuously links the superfluid and 
normal phases. One can easily show that hm = fJ-SF at 
H = H c i and fiM — fJ-N at H = H C 2, which ensures 
< y < 1 with y(H cl ) = 1 and y(H c2 ) = 0. 

For the discussions above, the mixed phase is assumed 
to be the ground state in the region H c i < H < H C 2- As 
a complete study, we have to prove that the mixed phase 
has the lowest energy in this region. Even though the full 
information of the scaling function Q{x) is still not clear, 
this can be done if the function Q{x) behaves sufficiently 
regularly. 

(1) The energy density of the mixed phase can be writ- 
ten as £ M {H) = ^£ f(H/H c i) with f(z) = z - |z 5 / 2 . 
From f'(z) = 1 — z 3 / 2 , £y[(H) is a monotonously de- 
creasing function of H in the region H c i < H < H C 2- 
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Combining Em = £sf at H = H c \ and the fact that £sf 
is iJ-independent, there is always £m{H) < £sf(H) for 
H c i < H < H C 2- 

(2) The condition £m(H) < £n(H) requires g(j) < 5(7') 

with g(t) = h/t - lg(t) (h/t) 5/2 , h = H/e F and 7' = 
H/hn(H). Even though we lack the full information of 
the scaling function Q (x) , it is sufficient to show £ m (H ) < 
£n(H) at H < H C 2, due to the continuity and the BCS 
instability £n(0) > £sf(0). From the first-order deriva- 
tive of g(t) at t = 7, g'(<y) = 7 " 2 ft [(H/H c2 f/ 2 - l], 
g(t) is a decreasing function near t = 7. Therefore, at 
< -ff c2 the condition 5(7) < g{"f') requires 7' < 7 or 
Mn > Mm- From ^n(O) > ^sf and A*m being an increas- 
ing function of H, the relation /ig-p < [1m < holds 
in the region H c \ < H < H C 2- Therefore, the condi- 
tion £m(H) < £n{H) is satisfied once the function G(x) 
behaves sufficiently regularly. A schematic plot of the 
energy densities for various phases is shown in Fig|2j 



yj 

UJ 




FIG. 2: A schematic plot of the energy densities (divided by 
the energy density of the ideal Fermi gas) for the superfluid, 
normal and mixed phases in the region H c \ < H < H C 2- 



III. RESULTS AND COMPARISON 



imbalance P c = P{H C 2) which can be determined in ex- 
periments and QMC calculations [28l. [29j. To this end, 
we consider a system with separately fixed and 
where the conversion between the two spin components 
is forbidden. From P(H cl ) = and P{H c2 ) = P c , the 
ground state is the unpolarized superfluid state at P = 
and the SF-N mixed phase for < P < P c . In the 
mixed phase, the effective "magnetic field" is given by 
H(P) = 7^€f(1 + 7P) 2 / 3 , and the critical population 
imbalance P r reads 



Pc = 



|S'(7) 



£-3/2 _ 2 g l(rf) 



(16) 



Note that this expression is consistent with the jump in 
Fig. [T] which can be calculated from Eq. (4). There- 
fore, the upper critical field H C 2 can be determined once 
P c is known. Further, the procedure can be simplified 
if we combine the expressions for the critical fields and 
the population imbalance, i.e. Eqs. (ITU)) and (fl"6|) . In 
this way, we obtain a very simple relation among the CC 
limits for different cases 



|^ = (1 + 7P C ) 2/3 . 



(17) 



For the case of fixed atom numbers, the energy density 
defined as £ = /i-^ri-j- + /i^n^ — V satisfies the relation 
£ = 3V/2 in all phases, since H is now no longer treated 
as an external field. To show the consistency between the 
theoretical formula of P c obtained above and the QMC 
calculations and cold atom experiments we derive the 
energy density £m as a function of the ratio n^jn^ [28| 



£ 



M 



. 3 (6^ 2 n t ) 2 / 3 /n,\ 



where the function I(z) can be shown to be I(z) = 
2-2/3^ [(i + 7 ) + (! _ 7 )z] 5 / 3 . The function I(z) is con- 
sistent with the formula used to obtain P c in the QMC 
calculations [28j . 

We can now determine the critical Zeeman fields H c \ 
and H C 2 from the known values of £, 7 and P c extracted 



Now we turn to determine the critical Zeeman fields 
H c i and H c2 . Since the assumption of fully polarized nor- 
mal phase is not correct for 7 < 5q and P c < 1, we may 
have to know the full information of the function Q{x) 
in the partially polarized normal phase. Fortunately, in 
the following we will find that H c \ and H C 2 can be com- 
pletely determined once the values of 7 and P c (as well 
as £) are known. 

From the exact expression (fTU[) . the lower critical field 
H c i can be determined once the values of £ and 7 are 
known. On the other hand, to determine the upper 
critical field H C 2, we need the value of ^'(7) which is 
not known so far. However, we find that we can de- 
termine the value of ^'(7) from the critical population 





7 Pc £ 


5'( 7 ) H cl [e F ] H c2 [e F ] 


QMC [28] 

MIT [291 
LKB-Lhomond [30] 
LKB-Lhomond [31J 


0.967 0.389 0.42 
0.95 0.36 0.42 
0.878 0.324 0.42 
0.897 0.359 0.41 


2.596 0.406 0.503 
2.464 0.399 0.485 
2.317 0.369 0.436 
2.586 0.368 0.443 



TABLE I: The data of 7, P c and £ from the QMC study [28j 
and experiments [29j - [3l| and the values of the critical Zeeman 
fields Hci and H C 2 determined from the formulae f| 10|) and 
(TTT)) . The value of G'ij) determined from Eq. |16]) is also 
shown. The value of P c in [3l| is not given, and we use our 
formula to extract it from the experimental data, see Fig. [3] 
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from QMC calculations and experimental measurements. 
The results are listed in Table HJ In Fig. [3] we show the 
data of the polarization P as a function of the Zeeman 
splitting H/ep extracted from the LKB-Lhomond exper- 
iments [31]. The critical Zeeman fields H c \ and H C 2 are 
around OAep and consistent with our calculations in Ta- 
ble U from the LKB-Lhomond data. The deviation be- 
tween the results from LKB-Lhomond data and those 
from QMC and MIT data comes mainly from the dif- 
ference in the values of 7. In the early studies [H, [29j . 
the value of 7 was reported to be around 0.96. However, 
recent data from LKB-Lhomond experiments show that 
this value becomes smaller, around 0.89. 




FIG. 3: (color online) The spin polarization P = (Nf — 
Ni)/(Nf + Ni) as a function of the external Zeeman field 
H (divided by ep) for a homogeneous Fermi gas at unitary. 
The blue dots are data extracted from the LKB-Lhomond 
experiments [3l|, |32| ■ The red solid line is the theoretical pre- 
diction. We use the data 7 = 0.897, £ = 0.41 from [HJ and 
H c i = 0.368eF. For the normal phase H > H C 2, we use the 
linear fit P = §xif/eF with \ = 0-54 suggested in [3]]]. In 



between, the formula P 



7 



\3/2 



1] is adopted. 



The value of P c is self-consistently determined as P c — 0.359. 



It was shown in 30, 31] that the normal phase of the 





Pc H c2 [e P ] 


X/Xo H [e F ] 


QMC [28] 

MIT [29] 
LKB-Lhomond [30J 
LKB-Lhomond [31] 


0.389 0.503 
0.36 0.485 
0.324 0.436 
0.359 0.443 


0.52 1.29 
0.49 1.35 
0.50 1.35 
0.54 1.23 



TABLE II: The spin susceptibility x (divided by xo) and the 
Zeeman splitting Ho at which the Fermi gas becomes fully 
polarized determined from different QMC and experimental 
data. 



unitary Fermi gas is well described by Landau's Fermi 
liquid theory. Here we can extract the spin susceptibility 
X of the normal phase from the determined upper critical 
field H C 2- The spin susceptibility x is defined by the 
linear relation between the spin polarization P and the 
Zeeman splitting H, i.e., P = (3xH) / (2xo £ f) where xo = 
3n/2ep is the spin susceptibility of an ideal Fermi gas. 
Therefore, the spin susceptibility \ can be determined 
by 



_X_ 

Xo 



Pr 



3 H c2 
2 e F 



(19) 



We can also determine the saturation splitting H at 
which the Fermi gas becomes fully polarized by setting 
P = 1. We find Ho/ep = 2xo/3x- The numerical results 
for x and Hq from different QMC and experimental data 
are summarized in Table HH We find good agreement 
among the results. The spin susceptibility of this uni- 
tary Fermi liquid is about half of that of the ideal Fermi 
gas, and the saturation splitting Hq is around 1.3ef- 

Finally, we compare our results with those obtained 
from the mean-field [33[ and beyond-mean-field [43[ the- 
ories. Our formulae (fTU| and (TTBl for the CC limits are 
model independent. In the mean- field theory, the uni- 
versal constant £ reads £mf = 0.5906 and the scaling 
function Q[x) takes the non- interacting form: Gmf{%) = 
\ [(1 + xfl 2 Q(\ + x) + (1 - x) 5 / 2 e(l - x)] . Numerical 
solution of equation © leads to 7mf = 0.8071 and there- 
fore <?mf(7mf) = 2.9307. Substituting these values into 
dlOl and fill we obtain 



H^ F = 0.693e F , P C MF = 0.933, (20) 



ijM F = 0.477e F , 



which agree well with the numerical values obtained 
in [33| . One finds that the mean field results of the criti- 
cal Zeeman fields H c \ and H C 2 deviate significantly from 
our results listed in Table HI The beyond-mean-theories 
may properly include the pairing fluctuations in the su- 
pcrfluid phase and the interactions in the normal phase 
and give reasonable values of the critical Zeeman fields. 
Such a calculation for the upper critical field H C 2 within 



the large-N expansion method has been performed in 43 1 
To the leading order of 1/N, the result reads 



H C 2 

£f 



0.693 



0.087 

N 



O 



1 

~N2 



(21) 



It is clear that the leading order result from the large-N 
expansion method is not consistent with our result, and 
the higher order calculations are needed. 



IV. INDICATION TO DILUTE NEUTRON 
MATTER 

For neutron matter the effective range of the nu- 
clear force, r ~ 2.7fm, is much smaller than the s- 
wave neutron-neutron scattering length, a nn ~ — 18.5fm. 
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Therefore, for dilute neutron matter which may exist in 
the crust of neutron stars, /cf^o can be relatively small 
but fcFdnn remains large. So the properties of the di- 
lute neutron matter are close to the unitary Fermi gas 
discussed in the paper. 

As a naive application, we can estimate the critical 
magnetic fields for the dilute neutron matter at which 
the superfluid state is destroyed and the matter becomes 
spin-polarized. The Zeeman splitting Ez = 2H in this 
case is H = fi n B where u n is the magnetic moment of 
the neutrons and B the magnetic field. After some simple 
algebras, we obtain the following formulae for the lower 
and upper magnetic fields 



n 



2/3 



-So, 



B r . 



7£(l + 7Pc) 2/3 (-) 



2/3 



B Q , (22) 



with no = 0.1 6fm 3 being the nuclear saturation den- 
sity and B Q = (37r 2 n ) 2/3 /(2M„/x„) ~ 10 19 G. In Table 
nH we calculate the critical magnetic fields B c \ and B C 2 
for some typical densities of dilute neutron matter from 
n/n,Q = 10~ 3 to n/no = 10 . We find that the critical 
magnetic fields are roughly in the range 10 16 — 10 18 Gauss. 
Therefore, the problem of imbalanced pairing and spin- 
polarization in dilute neutron matter is in principle rele- 
vant to compact objects known as magnetar s 1441 . which 
have surface magnetic fields of 10 14 — 10 15 G [45j . In fact, 
according to the scalar virial theorem which is based on 
Newtonian gravity, the magnetic field strength is allowed 
by values up to 10 18 G in the interior of a magnetar gi- 
ll. 



V. EXTENSION TO NONZERO 
TEMPERATURE 



The above model-independent approach can be gener- 
alized to finite temperature T, where both the normal 



neutron density 


n = 10" 3 n 


n = 10 2 no 


n = 10 1 no 


(frflnn)" 1 


-0.32 


-0.15 


-0.07 


Bd [Gauss] 


3.6 x 10 16 


1.7 x 10 17 


7.7 x 10 17 


B C 2 [Gauss] 


4.3 x 10 16 


2.0 x 10 17 


9.3 x 10 17 



TABLE III: The lower and upper critical magnetic fields, 
B c i and B C 2, for dilute neutron matter at different densi- 



ties, n/no = 10 



,10" 



,10 I . The values of (fea n 



also shown, from which we find that the effective couplings 
at these densities are really close to the unitary limit. In 
the calculations, we adopt the data £ = 0.41, 7 = 0.897 and 
P c ~ 0.36 from [3l[. 



and superfluid phases are spin-polarized due to the ther- 
mal excitations quasi-particles. From the universality, 
the EOS for the normal and superfluid phases read [§] 



Visi,sf(T,li,H) = 'Po(m)^n,sf ( — ,— 



(23) 



where we have set the Boltzmann constant fee = 1. The 
scaling functions for the normal and superfluid phases 
should be different. 

In the grand canonical ensemble, one expects that the 
phase transition along the T/ll axis is of second order 
at small H/fi and first order at large H/li. The first 
order phase transition is determined by the equation 
gn(H/fjL,T/n) = g SF (H/Li,T/Li), or explicitly H/li = 
W{T/ll) with known W(0) = 7. The first order phase 
transition should end at a so-called tricritical point 
(H/ fj,,T/fj,) = (a,b). At the mean- field level, it is pre- 
dicted to be (a, b) = (0.70, 0.38) H3- 

At fixed total particle number, fi is not a free 
variable, and the tricritical point is characterized by 
(Ttcp, -Htcp)- Due to the continuity with the zero tem- 
perature case, for T < Ttcp , there exist two critical fields 
H c i(T) = niW(T/m) and H c2 (T) = li 2 W{T/li 2 ), where 
Hi and fi2 are the chemical potentials corresponding to 
the superfluid phase at H = H c \ and the normal phase 
at H = H c2 , respectively. The region H c i < H < H c2 
for the mixed phase should decrease with increasing 
T, and finally disappear at the tricritical point where 
Htcp = aXrcp/fr and lli = li 2 = Trcp/b. 

When iV T and 2V; are fixed, for T < Ttcp, the SF-N 
mixed phase should be the ground state in the region 
Pi < P < P 2 with Pi = P(H c i) and P 2 = P(H c2 ). At 
T 5^ 0, Pi should be nonzero and increase with the tem- 
perature. At the tricritical point, Pi — P 2 — Ptcp is 
satisfied. Once the scaling function Q and the tricritical 
point (a, b) in the grand canonical ensemble are known, 
Ptcp and Ttcp can be calculated from the following 
model-independent formulae 



-Ptcp 
Ttcp 



\Q{a,b) - aG' x (a,b) 



bG'(a,b) 



= b 



5Ptcp 



2/3 



(24) 



with the definition Q' x [x, y) = dG(x, y)/dx an&Q' y {x,y) = 
dQ(x, y)/dy, where Q can be the scaling function of either 
the superfluid or the normal phase. 



VI. SUMMARY 

In summary, we have determined the lower and upper 
critical Zeeman fields for a homogeneous Fermi superfluid 
at infinite scattering length. Using the recent experimen- 
tal data from LKB-Lhomond, we found H c i ~ 0.37ep 
and H c2 ~ 0.44cf- The value of the lower critical field 
also gives a lower bound for the excitation gap Aq for 
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the balanced case. The results are highly related to the 
properties of dilute neutron matter in presence of a strong 
magnetic field which may exists in compact objects. The- 
oretically, it is interesting that we can prove the exis- 
tence of the two critical fields (H c i and H C 2) based only 
on the universal equations of state. We also presented 
a very simple proof for the fact that the mixed phase 
has the lowest energy, in contrast to the proof for the 
weak coupling case [20[. The value of the universal con- 
stant £ and the properties of the polarized normal phase 
are very important in calculating the CC limits. Fu- 
ture theoretical studies may focus on the calculation of 



the scaling function Q (x) from beyond- mean-field theo- 
ries [H, [43, Hil ■ Once the universal constant £ and the 
function Q (x) are known, one can directly obtain the crit- 
ical polarization P c from our model-independent formula 
(fT6|) and check the consistency between theories and ex- 
periments or QMC calculations. 
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